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Abstract 

The field-effect mobility of hole polarons in a-sexithiophene, measured in 
thin film transistors, was shown to be well fitted by Holstein's small polaron 
theory. Unfortunately, Holstein's formulation is based on an integral that does 
not converge. We show that the data are well fitted by a theory of polaron 
transport that was successful in accounting for mobility in molecular crystals 
of naphthalene. 



The prospect of field-effect transistors, FETs, made of easily processible thin organic 
films has stimulated considerable research effort. The most successful devices to date 
have been made of oligomeric rather than polymeric material, an outstanding one being 
a-sexithiophene, a-6T.0H From the device theory and operation it is possible to deduce a 
field-effect mobility, (Ifet, for holes moving in the active part of the device, the a-6T layer 
just above the gate dielectric. Recent measurements of mobility made in the space-charge 
limited current regime have been in good agreement with \ifet obtained in a long channel 
device,0 indicating that the device measurements can be considered to represent the mobility 
reasonably well. For convenience of notation we will drop the subscript in referring to the 
device measurements. 

The a-6T films for the FETs are usually made by vacuum evaporation onto the gate 
dielectric, the latter frequently being SiC>2 or MgF 2 . This results in a polycrystalline film 
which is, however, well oriented in that the long axis of the a-6T molecules stands almost 
perpendicular to the subtract S@ The transport direction in the FETs, parallel to the sub- 
strate, is then essentially perpendicular to the molecular chains. Mobility as a function 
of temperature has been measured in a-6T FETs by Horowitz et a/.,i Waragai et aM and 
Torsi et a/.@ In the two former cases the measurements were taken down to 100 and 77 K, 
respectively. Over this range mobility was seen to increase with increasing temperature. 
Waragai et al. attributed the transport to thermally activated hopping of polarons among 
the thiophene molecules. (Actually they used a derivative of a-GT, dimethyl sexithiophene.) 
Horowitz et al. analyzed the data down to 150 K in terms of free carriers, specifically holes in 
the valence band, undergoing multiple thermal trapping and release involving shallow traps. 
Below 150 K they attributed the mobility to the free carriers hopping among deep traps. 
Such a model could never account for the increase in mobility with decreasing temperature 
found by Torsi et al. below ~ 50 K.0 (See Fig. 1.) The model would be appropriate for 
carriers moving in wide bands of the kind found in three-dimensional semiconductors such as 
silicon, but not for carriers moving in the narrow bands characteristic of a molecular crystal, 
or the polycrystalline array of molecular crystals that a-6T samples generally consist of. For 
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these cases it is essential to include polaronic effects. 

Torsi et al. showed that their data are well fitted by Holstein's small polaron theory 
based on his Molecular Crystal Model.i In this model the electron (or hole) wavefunction 
is limited to a single site. The electron or hole becomes a small polaron, still limited to 
a single site, due to its coupling to the vibrations of a diatomic molecule at the site. It 
is well known that an excess electron or hole on a conducting polymer or on a conjugated 
molecule such as naphthalene or a-6T forms a polaron that is spread over a considerable 
number of sites, ~ 20 if the conjugation length allows. Although Holstein's calculations are 
readily generalized to such a polaron, and to the vibrations being the usual acoustic and 
optical modes of a lattice,! they are fatally flawed by the fact that their essential result, the 
integral obtained for the probability of a polaron moving or hopping from site to site, does 
not converge.00 Holstein obtains finite results for the high temperature (hopping) region 
and the low temperature (coherent motion) region by approximating the integral for each 
case with a leading term. Nevertheless, many of the ideas developed by Holstein and some 
predecessors form the basis for later theories of polaron transport. 

A form of the theory that was successful in describing mobility over a wide range of tem- 
peratures for naphthalene crystals starts from a simple Hamiltonian used by many workers 
in this field: 



where a' m creates an electron with energy e m at a site m in the crystal of N sites, 6 q creates a 
phonon with wave vector q and frequency co> q , V m - n is the intersite matrix element between 
sites m and n, g q is the dimensionless electron-phonon coupling constant and R m is the 
lattice vector that locates site m. The Hamiltonian is transformed to the polaron picture 
by a unitary transformation,!^ resulting in 
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+ £ V n _ m ele m A{A m + N- 1 ]T Ig^hu^ exp[iq ■ (R„ - K m )}AiAlA n A m , (2) 



n,m,q 



9 m = exp 



(3) 



q 

where v4 m creates a polaron at site m and 5^ creates a phonon of wave vector q that is 
associated with the motion of the ions about their displaced equilibrium positions. 

The Hamiltonian (2), (3) was used in the calculation of polaron mobility in naphthalene 
by Kenkre at al^ The formulation of KADD includes disorder, the site energy e m being 
taken as a randomly fluctuating quantity. The fluctuations can arise from crystal defects, 
which are surely present in these materials. It is usual to represent e m by a Gaussian or 
Lorentzian distribution of width ha. The latter was chosen by KADD. It was specified that 
a must be small compared to c<j q ; if a were not small a different approach to the problem 
would be necessary. 

As was done by Holstein and many others, KADD assumed that the hole interacts with 
only one phonon branch of mean energy Kujq and width hAi in the zth direction. The 
significance of this assumption will be discussed below. With these simplifications KADD 
obtained the mobility for field and current along the ith principal axis direction: 

flu = 2(ea 2 i )(k B Ty 1 \V t /h\ 2 ex P [-2g*coth(^(3huJo)] 

poo I 

x / rfte^%{2^sch(-/3^o)|^o(A 1 t)Jo(A 2 t)Jo(A 3 t)|[l + J 1 2 (A,t)/J 2 (A J t)] 1 / 2 }. (4) 
Jo 2 

where (3 = 1/ksT, the are dimensionless electron-phonon coupling constants as defined in 
Eq. (1), the <2j are lattice constants and the V{ the nearest-neighbor transfer integrals along 
the ith direction. The bandwidth W is taken as usual as 4Vi. It is seen that fin is similar to 
the result of Holstein's theory, but it converges due to the factor e~ at , representing disorder, 
in the integral. 

It is clear that many phonon branches, not only one, have significant interaction with 
the current carrier, whether electron or hole. For example, when a carrier jumps onto or off 
a chain, the well known distortion that occurs involves a number of phonon frequencies and 
more than one phonon branch. This was calculated explicitly for the case of polyacetylene.Q 



A similar polaronic effect must occur when a hole jumps onto or off an a-6T molecule. Of 
course the distortions are not large, but this represents a necessary energy expenditure that 
affects the hopping rate. A much bigger effect, however, can be caused by phonons that 
change the overlap of the molecules between which the polarons jump. For example, consider 
nearest neighbor hopping along that direction in an a-6T crystallite perpendicular to the 
long axis of the molecules where nearest neighbors are tilted at opposite angles to the normal. 
Librations around an axis in this direction would greatly affect overlap of adjacent molecules 
and therefore the rate of polaron motion. This appears to be the situation in at least one 
of the structures attributed to a-6T.0 In such a case the libration would have by far the 
strongest effect on the hopping and the effect of the other phonons could be considered to 
be included in the value of g^, which is determined by comparison with experiment in any 
case. 

As has been shown by many authors, in the approximation that only one phonon branch 
interacts with the carriers the thermally averaged bandwidth for the ith direction is0 

W T = A\Vi\ exp[-^ 2 coth(/?^ /2)] = 4|Vi| exp[-^ 2 (m + 1/2)] , (5) 

where m = [exp(/3h(Jo) — the number of phonons in the mode with frequency ujq. It is 
apparent that Wt decreases monotonically as the temperature increases from zero, rapidly 
when fcgT '/Tlujq > 1/2. Thus the mobility, which is proportional to W^, decreases with 
temperature due to this factor, rapidly at high temperatures. A second factor depending on 
T is csch(j3hu)o/2) in the integrand. This factor is proportional to m, the phonon abundance. 
In the low temperature limit this factor makes the argument of Jo small and therefore the 
integral small. As T increases this factor increases and the integral increases, quite rapidly for 
ksT > huj . Crudely it can be thought that the integral represents the hopping probability. 
Thus as T increases from zero the librons play the dual role of decreasing the coherent or 
bandlike motion through the decrease in bandwidth and of increasing the hopping rate. How 
strong these effects are depends on gi. 

Before comparing these predictions with the data of Torsi et al. shown in Fig. 1, we note 
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that their mobility values do not show the effect of trapping; their measurements have been 
made for gate voltages past the threshold value required to obtain a carrier concentration 
large enough to fill the traps. Also they have shown that grain boundaries do not dominate 
the transport, at least above the minimum in \i at ~ 45 K. Increase in grain size by an order 
of magnitude did not improve the room temperature mobility.0 From these facts Torsi et 
al. conclude that their FET mobility data represent intrinsic behavior of a-6T. Comparing 
the data of Fig. 1 with the theory discussed above, we see above the temperature of the 
minimum, Tl, a region in which /i increases rapidly with increasing T, consistent with 
hopping. Below Tl there is a region in which /i decreases rapidly with increasing T, as 
expected for the coherent motion. The fact that \x does not increase without limit with 
decreasing T, as predicted by Eq. ([|), must be due to another process, perhaps boundary 
scattering, determining \x in the low temperature limit. 

To fit the data of Fig. 1 with Eq. (f|) we took a% = 0.38 nm, the closest distance between 
adjacent molecules along the transport direction!! The widths of the phonon band, not 
known for a-6T, were taken from a libration of similar frequency which determines the 
polaron transport in naphthalene ,0 specifically hAi = 0.70 meV, HA 2 = 1.40 meV and 
hA 3 = 0.24 meV. In any case fin is not sensitive to the Aj values. The quantities \Vi\, gi, 
huo and a were treated as parameters. The fit shown in Fig. 1 was obtained for \Vi\ = 14 
meV, gf = 11.5, huj = 11.2 meV and a = 2 x 10 -5 ps _1 . The first three parameters are 
quite close to those obtained by Torsi et al. using Holstein' s results. 

It is of interest to compare these results with those obtained by KADD in fitting Eq. 
(f|) to the data for naphthalene. From 30 K, the lowest temperature for which there are 
data, to ~ 100 K /i for naphthalene decreases fairly rapidly with increasing T, although 
not as rapidly as for a-GT. Above 100 K fi continues to decrease with increasing T but the 
increase is much less rapid in two of the three directions shown. These results are consistent 
with Eq. (f|) with a somewhat larger huo but more importantly a smaller value of gi as 
compared with a-6T. Decreasing the electron-phonon coupling decreases both the rate of 
bandwidth decrease and the rate at which hopping grows with increasing T. The value of 
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gf for naphthalene, varying from 2.6 to 3.5 over the three directions, is small enough that 
the increase in the integral is not large enough to overcome the decrease in bandwidth and 
create a positive slope for mobility vs. temperature. Another contrasting feature of the 
naphthalene data is that the changes in slope of fi vs. T are much more gradual. This 
is apparently a result of a, the disorder parameter, being larger in naphthalene, 0.259 vs 
2 x 10~ 5 ps _1 . In fact, the later value for a is surprisingly small. We speculate that this 
may be a result of the holes being confined to within 1 or 2 molecular layers of the interface!! 
so that only a very small part of the a-6T layer is sampled, that part being more ordered 
perhaps because of the proximity of the interface. Of course, it is not possible to be certain 
of the parameters because there may still be some differences between [jlfet and the drift 
mobility, particularly below 45 K. 

If higher mobility is the requirement for improved organic FET performance, it is clear 
from Eq. (4) that for molecular crystals desirable properties are lower electron-phonon 
coupling and lower frequency of the interacting phonons. There is a limit to how far one 
can go in this direction, however, because polaron theory is no longer valid if g 2 hu>o/W < 1. 
From what was said earlier, naphthalene is superior to a-6T, its room temperature mobility 
ranging from 0.3 to 0.8 cm 2 /Vs in the different crystallographic direct ions. 

We acknowledge the support of the National Science Foundation under Science and 
Technology Center grant CHE912001. 
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FIGURES 

FIG. 1. hfet vs. T. Squares give experimental data of Torsi et al., 7 solid curve is the theo- 
retical fit described in the text. 
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